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We present the results of detailed dielectric investigations of the relaxation dynamics in 
DyMnOs multiferroic manganite. Strong low-frequency relaxation process near the paraelectric- 
ferroelectric phase transition is observed. We provide an experimental evidence that this relaxation 
mode corresponds to a chirality switching of the spin cycloid in DyMnOa. We demonstrate that 
the relaxation dynamics in DyMnOa is typical for an order-disorder phase transition and may be 
understood within a simple model with a double well potential. These results suggest the interpre¬ 
tation of the paraelectric sinusoidal phase in manganites as a dynamical equilibrium of magnetic 
cycloids with opposite chiralities. 


I. INTRODUCTION 

Multiferroic materials with a coupling of electric and magnetic degrees of freedom have attracted considerable 
interest after the discovery of a large magnetoelectric effect (ME) in several compoundsi^— . They are currently the 
subject of intensive study due to fascinating physical properties and potential for applications as multifunctional 
devices^i^. The rare-earth RMnOa manganites (R=Gd, Tb, Dy, Eu/Y) with orthorhombically distorted perovskite 
structure have emerged as a new class of multiferroics with strongly coupled antiferromagnetic and ferroelectric 
properties'*^. Several rare earth manganites orders antiferromagnetically below Tjy ^ 40 K into collinear paraelectric 
phased with sinusoidal modulations of spins. This phase is followed by a cycloidal spin order with nonzero electric 
polarization below Tq ^ 20 K. 

It should be noted that the Mn^+ atom is supposed!^ to have a Heisenberg spin with a fixed length of S' = 2. The 
purely sinusoidally modulated spin phase contradicts this property. Model calculationsiS have obtained the sinusoidal 
order only as a time-space average of the simulated cluster. Thus, there is a possibility that a short-range dynamic 
order exists in the intermediate temperature range Tc < T < Tjv which is responsible for the “hidden” spin. As 
the magnetic order at low temperatures is spin cycloid, it is natural to assume that the dynamic short-range order is 
also a spin cycloid. This would imply that there are fluctuating ferroelectric regions in the sinusoidal phase and that 
the ferroelectric transition is actually of the order-disorder type. Such transition has been also suggested in Ref. m 
where a c-axis relaxation typical for the order-disorder type transitions has been investigated. The fact that the wave 
vector qMn of the spin wave does not change at the transition temperature Tc [ll] is also an indirect evidence that 
the ferroelectric transition is not of the displacive type. 

The data by terahertz spectroscop y^^d^ evidence the nonzero dielectric contribution of electromagnon in the sinu¬ 
soidal phase. According to the commonly accepted mechanism of the electromagno n^ ^ the majority of the spectral 
weight of this mode originates from exchange striction mechanism and can only exist in magnetic phases with non- 
collinear spin arrangement. These facts again favor the hypothesis of dynamical cycloidal spin order in the sinusoidal 
phase. 

Recent theoretical analysis of the terahertz dynamics in the sinusoidal phase suggested an explanation based 
on anomalous magnetoelectric coupling. Investigations of the collinear sinusoidal phase in the diluted compounds 
TbMni_ 2 ,Ab 03 [13 and the observation of the memory effect in the low temperature sinusoidal phase in the multifer¬ 
roic MnW 04 [13 have suggested the presence of the nanosize ferroelectric domains and support relaxor order-disorder 
type transition. 

In this work we present the analysis of the critical behavior of the low-frequency relaxation in DyMnOa. The 
observed critical behavior confirms that the sinusoidal to cycloidal phase transition is of the order-disorder type. Our 
model suggests the presence of the short-range cycloidal order in the collinear spin phase. 


A. DyMnOs 

In zero magnetic field, DyMnOs undergoes an antiferromagnetic (AFM) transition with a temperature dependent 
modulation vector (0 qMn 1) around Tjsf = 39 K. With further cooling, the qMn value is locked at the transition 
temperature, Tc ~ 19 K, and ferroelectric (FE) polarization appears simultaneously along the c-axis. It is well 
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established that the ferroelectric polarization in DyMnOa is induced by a cycloidal magnetic orde r^d°d^ through the 
inverse Dzyaloshinskii-Moriya (DM) interaction^i^^ri^i and it can be written as 


P = Y,Kei 


X {S, X S,) 


( 1 ) 


where eij denotes the unit vector connecting the spins Si and K ]s & constant representing the exchange interaction 
and the spin-orbit interaction. Accordingly, the electric polarization, P, is intimately linked to the chirality of the 
magnetic cycloid (clockwise respectively counterclockwise cycloidal magnetic ordering), that is, changing the direction 
of +P —?> —P implies changing the rotation (chirality) of the magnetic cycloid. This was demonstrated^^ by the 
asymmetry in the scattering of left-hand and right-hand circularly polarized x-rays by nonresonant magnetic x-ray 
diffraction for closely similar compound TbMnOa. 

The unique cycloidal magnetic ordering below Tc in RMnOs perovsikes is assigned to the Mn-3d spinsii^i^^— but 
ordering of the Dy-4f moments is also of interesli^ '^^i^^ . Although it is the magnetic structure of Mn subsystem that 
determines the emergence of ferroelectricity in rare-earth manganites^S, the mutual coupling of the Mn-3d and Dy-4f 
moments and, consequently, the ordering of the 4f moments in Dy causes a particular large polarization observed in 
this materiali^i^. In addition, the basal plane of the spin cycloid flops from the be plane to the ab plane and rotates 
the polarization by applying a magnetic field along the 6-axis. 

Concerning the temperature range Tc < T < Tjv, there have been a number of studies of the magnetic structure 
of DyMnOa and the related material TbMnOs using techniques such as neutron diffraction^i^^ and magnetic x-ray 
scattering2ii^. These studies have all shown that there is a long-range sinusoidal magnetic ordering of the Mn^+ ions in 
this temperature range, with the moments aligned parallel to the crystallographic 6—axis, although Mannix et al. 
and Wilkins et al. [33|, both show a small component of the moment is aligned with the c—axis. 


24| 


II. EXPERIMENT 

DyMnOa single crystals were grown in Ar flow by a floating-zone method with radiation heating^^*^. Terahertz 
properties of the samples from the same batch have been presented prevlously^S*^. The complex dielectric constant 
was measured for electric field along the crystallographic a-axis in the frequency range O.lHz-lMHz using a frequency 
response analyzer in magnetic fields 0-14 T and with i3||6-axis. We investigated the temperatures near the critical 
temperature of the paraelectric-ferroelectrlc phase transition (red circles in Fig. [5]) with an increment of AT = 0.1 K 
by using a Physical Property Measurement System (PPMS). Silver paint contacts were applied to the sample forming 
a capacitor. 


III. PSEUDO-SPIN-MODEL 


The theory of dynamic critical phenomena for a multivariable system can be formulated as a generalization of the 
single particle Langevln equation. The equation of motion for the time-dependent local configuration of the order 
parameter field is most conveniently given by the time-dependent Ginzburg-Landau equation^^— . 


dP[r,t) 


6F 


SP(r, t) 


C(i’,0 ■ 


( 2 ) 


Here, the order parameter is given by static electric polarization P, P is a dissipation parameter, and is a 

random noise simulating the effect of thermal excitation of the order parameter. In order to guarantee that the 
system reaches the canonical equilibrium probability distribution at long times, C,{Y,t) is a random Gaussian variable 
satisfying = 0 and = 2DkBT5(v — v)5(t — t) [^, where D Is the diffusion coefficient. The 

crucial point is to find an expression for the free energy in Eq. ([5]) . In the framework of the Landau theory for a second 
order phase transition, Smolenskii^ introduced a bi-quadratic term (P oc yP^M^, magnetodielectric effect) into the 
free energy which accounts for the coupling between magnetization and electric polarization. Bi-quadratic terms are 
invariant to all symmetry operations and thus they are allowed In any material with coupled spin and charge degrees 
of freedom. Since the dielectric susceptibility is determined by taking the second derivative of the free energy with 
resp ect to the polarization, the dielectric constant will be proportional to the square of the order parameter, e oc 
|43l| . Describing magnetodielectric effects in antiferromagnetic materials, the expression F cx yP^M^ is not sufficient 
since the magnetization, M, remains zero in the ordered phase. In such a case, M is replaced by the antiferromagnetic 
vector, L = Mi — M 2 . Here Mi and M 2 are the magnetizations of two antiferromagnetic subsystems. 
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Within a more general model, Lawe a'^^d^ et al. proposed the coupling of the polarization to the g-dependent 
magnetic correlation function {MgM_q). This coupling leads to a magnetodielectric term in the free energy F oc 
9iQ)P^{MqM-q), where g{q) is a g-dependent coupling constant. The g-dependence of the free energy via a spin- 
spin correlation function enables to apply it to very general forms of magnetic order, including ferromagnetic (FM) 
and antiferromagnetic (AFM) transitions. In order to obtain a microscopic theory for g(q) in systems with a strong 
spin-lattice interaction, the coupling between the polarization and the spin correlations arises from the coupling of 
magnetic fluctuations to the optical phonons. That is, the spin correlations perturb the optical phonon frequencies 
which in turn shift the dielectric constant through the spectral weight transfer and the Lyddane-Sachs-Teller relation. 
The model determines the coupling g(q) by expanding the exchange integral of neighboring spins in terms of the 
normal coordinates for the phonons. Physically, this procedure corresponds to a coupling between the magnetic 
correlation function and atomic displacements. 

In multiferroic rare earth manganites, RMnOa, the electric polarization is directly linked to the chirality of the 
magnetic cycloid^iS^. Based on this fact, we propose that the polarization in the ferroelectric phase in DyMnOa is 
proportional to the difference of opposite chiralities of Mn^+ magnetic cycloids. Here we assume an order-disorder type 
phase transition between paraelectric and ferroelectric states. Similar analysis in a triangular lattice antiferromagnet 
RbFe(Mo 04 ) demonstrated a proportionality between polarization in the multiferroic phase and the chirality difference 
of the magnetic structure^. 




t” 


1’’ T>Tc 



CfR ‘^R+l 



FIG. 1: Pseudo-Spin Model (a) Short range cycloidal order of the Mn®'*' magnetic moments for T > Tc- (b) According to 
Eq. 0, electric dipoles associated with the ions are generated by the canting of neighboring spins leading to a mesoscopic 
electric polarization (red and green zones), however, the macroscopic electric polarization is zero for T > Tc- Each electric dipole 
interacts with neighboring electric dipoles, (c) Ising type pseudo spins in local double-well potential separating energetically 
the clockwise a6-cycloidal and the counterclockwise a6-cycloidal order of the Mn®"*" magnetic moments. Each pseudo spin is 
in a local double-well potential and interacts with neighboring pseudo spins by harmonic forces, represented as springs. Each 
pseudo spin generates a electric dipole moment of pR = auqAx. Although the microscopic mechanism of the magnetoelectric 
coupling is more complex, within present simple model it is represented as an interaction between neighboring pseudo spins 
via harmonic forces, (d) Long range cycloidal order of the Mn®"*" magnetic moments for T <C Tc and leading to (e) non-zero 
macroscopic polarization, (f) At low temperatures most of the pseudo spins occupy the same side of the double-well potential. 


In the present model the following assumptions for the phase transition in DyMnOa are imposed: (i) A disorder 
between clockwise and counterclockwise a6-cycloidal order of the Mn^+ magnetic moments is assumed, (ii) The electric 
dipole moments are associated with the displacement of the ions due to inverse DM interaction^i^S, (iii) The 
direction of the electric dipoles depend on the chirality of the magnetic order, (iv) Two possible direction of the 
electric dipoles are energetically separated by an energy barrier, (v) Similar to Lawes^^f^ et al. we propose a coupling 
of the magnetic correlation function and the correlation of atomic displacements. Thus we assume that we can 
describe the ordering of the magnetic sublattice by the ordering process of the ions. 
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A useful tool for studying phase transitions is provided by the unidimensional single ion model^^— . This model 
contains an array of atoms linked by harmonic forces, with one atom in each unit cell (Figll}. Each atom is located 
in its local double-well potential which represents the rest of the crystal. 

Assuming order-disorder limit ^ 1 the system can be described by the pseudo spin formalism 

and thus the model Hamiltonian is essentially governed by an Ising type interaction in combination with an interaction 
of dipoles with a homogeneous electric field^^— 

H (X — E Aa: JR'^R ^ ^ E{t)q/S.X(jR , (3) 

R,R R 

where /S.x is the displacement of the ion, ^ is the coupling constant between ions at position R and R, 
E(t) is a time dependent homogeneous electric field, q is the charge of the oxygen ion, and an is the pseudo spin at 
position R with an = xr/A x which can take the values -fl and -1. 

Introducing a statistical mean value of the spin variable, s = {N^ — N_)/{N^ + N_), with A+ and N_ being the 
occupation numbers of ions in -|-1 und -1 state respectively, and neglecting fluctuations of correlations, the free 
energy of the system with polarization P = nqAxs is given by^ 


F = -ksTlnZ = / i - nksTln 

j I nq-^Ax 

where n is the dipole density and Jq = J2r Jr r represents the coupling constant characterizing the interaction of an 
0^~ ion at position R with another ion at position R located within an interaction radius. Since the first term 
in Eq. ([3|) is bilinear in the displacements of ions it corresponds to a electric dipole-dipole interaction. 

Expanding P with respect to P and E and taking the functional derivative of F with respect to P, leads to an 
equation of motion of the homogeneous order parameter: 


2 cosh 


^^P + E{t)qAx' 


dV 


( 4 ) 


dPjt) _ SF 
dt SP 


f Jo 2(Aa:)Vg \ 32 (Ax)V 4 ^3 2(Aa:)Vo 

nq’^ksT J Gn^q'^k^T^ ksT 


( 5 ) 


Here T = nq^J^ ^voexp{—E a/ kRT) is the damping parameter. 

According to Landau theory^ the phase transitions takes place when the term in Eq. linear in P vanishes. 
This determines the phase transition temperature as Tq = 2Jo{Ax)‘^/ks- Substituting E{t) = SEexp{—iujt) and 
P{t) = Ps + SP exp{—iojt) [s^, where Pr is the static electric polarization, and comparing the relaxation rate in 
a single double well potentia l^'^’^^i^^ we get an expression for the relaxation strength Asr = SP/{eoSE) and the 
relaxation time r in close agreement with Lines and Glass and Blinc and Zeks as 


Ae 


-1 

r 


Fi^iT-Tc) T>Tc 
T<Tc 


( 6 ) 


(27rT) ^ 



exp {-^) 

f (2^) exp (-^) 


T>Tc 
T <Tc 


( 7 ) 


Here vq is the attempt frequency and Ea is the energy barrier separating two local minima (Fig. [T]). 

Equations ([6]) and 0 correspond to an order parameter which is homogeneous throughout the entire volume of the 
crystal. For a spatially dependent order parameter, P(r), and assuming the Gaussian approximation, the relaxation 
time in Fourier-space become a^^'^^ . 


T{q) 1 = (2a -f bq'^) T , (8) 

where a = A(T — Tq), 6 is a constant and q is the wave vector. As a consequence, we see that each Fourier component 
of the order parameter behaves as an independent particle connected to a spring^S and the fluctuations in each mode 
decay with a different relaxation time. Only in the long wavelength limit, g = 0, the relaxation time in Eq. m is 
equal to that in Eq. ©• 
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In the opposite (displacive) limit of the model, EA/ksTc <C 1, the well-known expressions of Ginsburg-Landau 
mean-field theory are obtained for g —>■ 0: 


Ae 


-1 

r 


(^{T-Tc) T>Tc 
l^(Tc-T) T<Tc 


(9) 


(27rr) ^ 


^{T-Tc) T > Tc 
{Tc -T) T <Tc 


( 10 ) 


Here A = 30fcs/(87rJo), Tc = rj/A, and C are constants of the displacive limit; 0, and rj are microscopic parameters 
characterizing the nature of the potential in the cjA single ion model^^— . Nota bene, in canonical ferroelectrics the 
ferroelectric phase transition of the displacive type is accompanied by softening of a characteristic phonon. In the long 
wavelength limit an overdamped softening mode is characterized according to Eq. m and Eq. Recently, a 

critical slowing down at the ferroelectric phase transition has been observed in a chiral multiferroic MnW 04 [^. This 
behaviour has been attributed to a overdamped softening of an electromagnon mode^i^ obeying the temperature 
characteristic given in Eq. m and Eq. (nni with a critical exponent larger than 1. 


IV. RESULTS AND DISCUSSION 


As demonstrated in Eig. [2^, the dielectric permittivity closely follows the known phase diagram of DyMnOa for 
R||6-axisi. The changes in the dielectric permittivity are especially strong at the transition to the ferroelectric state 
with P\\a (Eig.Oa)- In addition a well pronounced absorption is observed in the close vicinity of Tc (Fig.j^t) Detailed 
analysis of the low-frequency dielectric relaxation is presented for magnetic fields of 4T,10T and 12T corresponding 
to the transition from the sinusoidal paraelectric to the P||a-axis ferroelectric phase. 

Figure [3] shows typical dielectric spectra of DyMnOs close to the ferroelectric transition temperature Tc ~ 18 K. 
The spectra below the megahertz range are dominated by two relaxation processes. Only a wing of the high-frequency 
relaxation is seen in our spectra because the characteristic frequency of this mode is far above 1 MHz. According to 
previous dielectric studies, the high-frequency mode can be attributed to the relaxation of the domain walls^^^ in 
DyMnOs. 

In the following we concentrate the analysis on an absorption peak observed for frequencies below 1 kHz. As 
seen already in the spectra in Fig. (imm. this peak grows in magnitude with decreasing temperature, reaches a 
maximum value at Tq, and decreases again after passing the critical temperature. The observed dielectric relaxation 
is slightly assymmetrical with broadening towards low-frequencies. 

To obtain quantitative information of the origin of the low-frequency mode the spectra were fitted to the phe¬ 
nomenological Havriliak-Negami equation^: 


e' -h e" = 


Ae 


(1 -I- {iojTy °‘)P 


( 11 ) 


Here Ae is the relaxation strength, Soo is the high frequency limit of the dielectric constant, r is the characteristic 
relaxation time, a and /3 are the width and asymmetry parameters, respectively. A simple Debye behaviour in Eq. (ED 
would corresponds to a = 0 and /3 = 1. Values of 0 < a < 1 and 0</3<(l — a) ^ results in broadened asymmetric 
loss peaks with power laws of and as the low- and high-frequency asymptotic behaviour, respectively. 

Since the symmetric Cole-Cole function given by 


e' -f e" = 


Ae 


(1 -I- (iwr)i “) ’ 


( 12 ) 


is intensively used as a fitting function to describe permittivity data of conventional materials as well as magnetoelectric 
materialsiii^, we demonstrated by the dashed line in Fig. [3] for the T = 16 K data, that the symmetric Cole-Cole 
function results in a worse fit to the data compared to the Havriliak-Negami expression. Therefore, the subsequent 
analysis within the present work has been done according to Eq. ED- Most probably, a fitting procedure using the 
Cole-Cole function would lead to qualitatively similar behaviour of the relaxation time and dielectric strength. 

The deviations from a Debye spectral shape of the relaxation are commonly ascribed to a distribution of relaxation 
times. According to Eq. (|S]), this can be caused by fluctuations of the order parameter near the phase transition 
temperature. The mean logarithmic relaxation time is related to the characteristic relaxation time by^ 


{InTHN) = Inr -I- 


ijj{P) + Eu 


1 — a 


(13) 
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FIG. 2: Phase diagram of magnetoelectric DyMnOs in magnetic fields along the b-axis. (a) The color map is a contour plot of 
the real part of the dielectric constant along the a-axis measured at 10 Hz. The ground state of DyMnOa in zero magnetic field is 
a 6c-cycloid with P||c. In external magnetic fields the electric polarization flops from P||c to P||a. (b) Low frequency dispersion 
phenomena near Tc- (c) Logarithmic plot of the imaginary part of the permittivity reveals a well pronounced absorption in the 
vicinity of the paraelectric to ferroelectric phase transition and the signature of domain wall relaxation phenomenon at higher 
frequencys^. 
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FIG. 3: Spectrum of the dielectric permittivity in a magnetic field of fj.oH = 12 T. Typical frequency dependence of the dielectric 
constant, Re(e), and dielectric loss, Ini(£), at several temperatures above (T = 18K, T = 20K) and below (T = lAK, T = 15R', 
T = 16K) the critical temperature, Tc « 17 K (at 12 T). Black solids lines represents the fit according to the Havriliak-Negami 
relaxation function, Eq. m- The dashed black line shows the ht according to the Cole-Cole equation, Eq. m- Dash-dotted 
lines with open circles schematically indicate the position of the spectra with increasing temperature. 


where ipiP) is the digamma function and Eu ~ 0.577 is the Euler constant. The width of a non-Debye relaxation is 
defined as the variance, cr^ = {(Iuthn)^) — (lnr//jv)^, of the distribution of logarithmic relaxation times and for a 
Havriliak-Negami function is given by^ 


,2 ^ V’'(/3) ^ 

(l-a)2 6(l-a)2 3 


(14) 


In present work we do not consider the relaxations by the domain walls. Compared to fluctuations on the atomic 
level, ferroelectric domains are typically^ of /rm size and they are responsible for the high-frequency dielectric relax- 
atioii^i^. Further two arguments are in favor of nanosize origin of the relaxation discussed here: (i) the low frequency 
relaxation (Fig.|31) is well pronounced below and above Tc^ and (ii) no signature of thermally activated creep motion of 
domain walls is evident in the Cole-Cole plots (Fig. [d]). In the latter case a linear relationship between the imaginary 
and real part of the permittivity with Im(e)oc tan(7r/?/2)-Re(e) and 0 < ,5 < 1 is expected^^— . 


A. Relaxation Strength 

The temperature dependence of the inverse relaxation strength of the low frequency mode is presented in Fig. [S] 
Dielectric permittivity diverges as the temperature approaches Tc- At the ferroelecric transition temperature 
DyMnOa undergoes a phase transition accompanied by a minimum in at Tq- As predicted by Eq. (p)) . 
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FIG. 4: Cole-Cole plots of low frequency relaxation. The inset shows the temperature dependence of the relaxation time in the 
Arrhenius representation (see section HV Bl) . 



12 14 16 18 20 22 
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FIG. 5: Relaxation strength in DyMnOa. Inverse relaxation strength along the a-axis of the low frequency mode as a function of 
temperature for different magnetic fields along the 6-axis. The yellow symbols are obtained by the spectral analysis procedure 
and the red solid line correspond to fit function according to Eq. ®. 
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FIG. 6: Crystal structure of the orthorhombic rare earth manganite DyMnOz. Green spheres are Dy ions, blue spheres represent 
Mn ions and magenta spheres are O ions in Pbnm setting. Solid gray lines mark the unit cell. Lower frame shows the ab-plane 
of DyMnOs cut at the Mn02 plane (dashed gray lines). 


close to transition Ae~^{T) shows the regions of linear dependence. This behaviour demonstrates critical dynamics 
characteristic for order-disorder phase transitions. 

A significant rounding near Tc is attributed to fluctuations of the order parameter and to limitation of the correlation 
length close to Tq [1^ . In addition, the straight lines of the model fits cross the x-axis at temperatures deviating 
by about one Kelvin. These effects are not accounted within the present simple model since it implies a molecular 
field approach and neglects fluctuations. Alternatively, a narrow temperature range near Tc may be analyzed using 
the formalism of critical exponents oc |T — Tcl"^■ Such analysis (not shown) reveals critical exponents for the 

relaxation strength of the low frequency mode near Tc between 7 = 1.8 and 7 = 2.3 depending on the magnetic field. 
These rather high values of the critical exponents may also be seen as an evidence for a order-disorder type phase 
transitions^^. 

The electric dipole moment Is generated by the displacements of the ions and only the Ions labeled as 02 
In Fig. ini generate a static electric polarization along the a-axls. Therefore, we assume an electric dip ole density in 
DyMnOa equal to 2/3 of the density of Ions in the unit, no 2 - = uq ^ 3.5 • [TIIIt^ and a charge 

of g = 2e, where e is the elementary charge. According to Eq. the slope of Ae^ ^ Is Inversely proportional to 
n{qAxy = 2Pg /n with static electric polarization along the a-axis Pg = nqAxj^/2. The factor \/2 appears to to 45° 
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degree misalinement between the oxygen displacement and the a-axis (Fig. [ 6 ]). Hence, for T < Tq from the slopes 
of we may directly estimate the static electric polarization. The microscopic parameters of the model, Ax 

and Jo are obtained from the values of the static polarization and from Jq = ksTc/2{Ax)'^. These parameters are 
given in Fig. [5] Taking into account the simplicity of the model, the obtained values of the electric polarization agree 
reasonably well with directly measured datai P ^ 2000^C/m^. In addition, the obtained coupling constant coincide 
with the classical electric dipole energy in the Mn-O-Mn chains along the b-axis Jo ^ {q^ /2Treo)/ {P + ^ 70J/m^ 

(Fig. ini with I = 2.22A and s = 1.90A being the long and short bond distances^^ . 

The temperature dependence of the inverse relaxation strength becomes steeper in the paraelectric phase T > Tq 
although the model predicts an opposite behavior. This effect is evidently not captured within the assumptions of 
present simple model. Two possible explanations for this behaviour may be suggested in present stage: (i) Change 
in effective dipole density at the phase transition. In the disordered state a large portion of the magnetic cycloid is 
distorted and is included into the border regions between the left- and right-rotating cycloids. The oxygen ions in 
these regions are effectively excluded from the relaxation process thus reducing the effective dipole density n of the 
model. This explanation has been used in the fits to the data in Fig. [5] and the effective dipole density for T > Tc 
is indicated at the fit lines, (ii) The model can be modified to account for higher order terms^^i^i. A coupling term 
between the polarization and magnetization may me explicitly included into the free energy expansion, Eq. ([S]). The 
term of which is always allowed by symmetry is Thus the free energy expansion near Tc, Eq. 0 , can be 

written as 

F = aP‘^ + bP^ - gEP + , (15) 

where a = A{T — Tq), and b, g, and 7 are the constants of the model. In these modification, M has to be understood 
as the amplitude of the transverse component of the spin-cycloidal S = {0,M cos qy, M sin qy). Applying Eq. (fT5l) to 
Eq. ([2]) leads to 


dP (t) 
dt 


-F {2aP + AbP^ - gE + 'yPM^) 


(16) 


With E = dEexp{—iojt) and P = Ps + 6Pexp{—iojt) and assuming that oc {Tc — T) = —Ca below Tc and 
= 0 above Tc the relaxation strength and relaxation time becomes 


Ti^iT-Tc) T>Tc 
^ 12(1- 7 ^) 7 ^^ (Tc-T) T<Tc 


(17) 


(27rr) ^ = 


I 

Vo (T 

TT V 

[ 2 ( 1 - 

yC) 


exp 


V keT J 


{ Ea 

V fcsT 




T>Tc 
T <Tc 


(18) 


It is clear, that adding explicitly magnetoelectric terms into the free energy, the ratio of the slopes of the relaxation 
strength below and above Tc deviate from the value of 2 as predicted by a simple model. In this case, the slope below 
Tc is governed by a new parameter (1 — yC) in Eqs. (I17I18I) . 


B. Relaxation Time 

The inverse relaxation time of the low frequency mode obtained from the spectral analysis is shown in Fig. [T] On 
cooling from the paraelectric phase the inverse relaxation time decreases toward Tc- Below Tc, 1/r shows a broad 
characteristic maximum and decreases again for low temperatures. Qualitatively, the temperature behavior of the 
relaxation time can be explained as a superposition of two processes: (i) Activated behaviour with a characteristic 
energy of Ea 100 K and (ii) Critical slowing down of the relaxation in the vicinity of Tc- This observation is typical 
for order-disorder phase transitions involving shallow double well potential s^^i^^F^ . 

Both processes determining the temperature evolution of the relaxation time are qualitatively well captured within 
the present simple model. The temperature activated behaviour, expressed by the exponential factor in Eq. 0 , 
prevails for temperatures far from Tc and, therefore, causes an overall decrease of the relaxation time for decreasing 
temperature. Qualitatively similar behaviour of the relaxation time was found^i for a c-axis relaxation (e||c) in 
DyMnOa for a transition to a 6 c-cycloidal magnetic ordering. The P||c-axis state is achieved in Dy for cooling in 
zero external magnetic field. It seem to be plausible that the c-axis relaxation in multiferroic manganites may be 
explained by the present model as well. 
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FIG. 7: Relaxation time in DyMnOs. Inverse mean relaxation time along the a-axis of the low frequency mode as a function 
of temperature for different magnetic fields 731| 6-axis. The yellow symbols are obtained by the spectral analysis procedure and 
the red solid line correspond to fit function according to Eq. 0. 


In case of inverse relaxation time the suggested model gives a qualitative explanation of the observed data. Eq. © 
contains two temperature-dependent factors, Arrhenius term exp{—EA/{kBT)) and the critical-slowing term voiT — 
Tc)/T. These two terms qualitatively explain the temperature dependence of the relaxation time close to phase 
transition. In order to obtain reasonable fits to the critical behavior of the relaxation time, different values of 
the attempt frequency above and below Eq have been used as well as a constant was added to Eq. ©. Nota 
bene, accounting magnetodielectric effects the slope of the inverse relaxation time is not necessarily 2, see Eq. m- 
Thus the values of the attempt frequency above and below Tc can take the same value. Eventually, in addition to 
magnetodielectric effects, the distribution of relaxation times plays also an role to explain that feature. 

The ratio of EA/ksTc ~ 5 is a further evidence of an order-disorder type phase transition^ and is in contrast to 
a displacive type phase transition where EA/ksTc <C 1 holds. Furthermore, EA/ks ^ 100 K corresponds well to a 
characteristic energy of the magnetic order as determined by Neel temperature of Tjy ^ 39K (s^ . 

Finally, Fig. [8] shows the width of the low-frequency dielectric relaxation in DyMnOa, obtained via Eq. (fl3 . The 
increase of the characteristic width below Tq ~ 18 K is clearly seen. In the ordered magnetic state the effective length 
of the elementary cycloids increases thus leading to broader length distribution. Most probably, this also leads to the 
observed broadening of the dielectric relaxation. 


V. DISCUSSION 

The analysis of the low-frequency dielectric relaxation demonstrates an overlap of two processes: critical dynamics 
close to Tq and activation behavior in the broader frequency range. This observation may be reasonably explained 
using a simple model of an order-disorder phase transition with a double well potential. This potential reflects a 
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FIG. 8: Relaxation width in DyMnO^. Asymmetry and width of the low-frequency dielectric relaxation in DyMnOa, obtained 
via Eq. (ITU 


dynamical switching between cycloids of the opposite chirality. Several parameters of the model correlate well with 
physical properties of DyMnOa. Thus, the characteristic energies of magnetic ordering and the value of the static 
electric polarization are in agreement with known values. We note that present experiments could be analyzed up to 
T = 22 K only. Strictly speaking, the dynamical behaviour of the cycloids still may change close to T^. 

Most importantly, the present experimental data and the simple model suggest to explain the paraelectric sinusoidal 
phase in rare-earth manganate as a dynamical equilibrium of cycloids with opposite chiralities. In addition to the 
dielectric results, this hypothesis resolve several experimental constraints which contradicted the concept of static 
sinusoidally modulated magnetic phase. 


VI. CONCLUSIONS 

Low-frequency relaxation mode is observed in dielectric properties of DyMnOa multiferroic manganite and it reveals 
critical behaviour at ferroelectric transition temperature, Tc 18 K. Together with temperature-activated relaxation 
rate, the observed mode may be qualitatively explained within a model for order-disorder phase transition. The 
model assume a switching between magnetic cycloids with opposite chirality and correlates well with known physical 
properties of DyMnOa. 

Combining present results with several other experiments on multiferroics we suggest that the paramagnetic sinu¬ 
soidal phase should be explained as a dynamical equilibrium between the clockwise and counterclockwise cycloidal 
magnetic orders. The short range order in the paraelectric phase is transformed to a long-range cycloid at the 
ferroelectric transition temperature. 
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